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Abstract 

We consider a factorization of the non-stationary Schrodinger ope- 
rator based on the parabolic Dirac operator introduced by Cerejeiras/ 
Kahler/ Sommen. Based on the fundamental solution for the parabolic 
Dirac operators, we shall construct appropriated Teodorescu and Cauchy- 
Bitsadze operators. Afterwards we will describe how to solve the nonlinear 
Schrodinger equation using Banach fixed point theorem. 
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1 Introduction 

Time evolution problems are of extreme importance in mathematical physics. 
However, there is still a need for special techniques to deal with these problems, 
specially when non-linearities are involved. 

For stationary problems, the theory developed by K. Giirlebeck and W. 
SproBig [TU], based on an orthogonal decomposition of the underlying function 
space in terms of the subspace of null-solutions of the corresponding Dirac op- 
erator, has been successfully applied to a wide range of equations, for instance 
Lame, Navier-Stokes, Maxwell or Schrodinger equations [5], [TU], [H], [2] or 
[13] . Unfortunately, there is no easy way to extend this theory directly to non- 
stationary problems. 

In [7] the authors proposed an alternative approach in terms of a Witt basis. 
This approach allowed a successful application of the already existent tech- 
niques of elliptic function theory (see [lOj . [5]) to non-stationary problems in 
time- varying domains. Namely, a suitable orthogonal decomposition for the 
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underlying function space was obtained in terms of the kernel of the parabolic 
Dirac operator and its range after application to a Sobolev space with zero 
boundary- values. 

In this paper we wish to apply this approach to study the existence and 
uniqueness of solutions of the non-stationary nonlinear Schrodinger equation. 

Initially, in section two, we will present some basic notions about complexi- 
fied Clifford algebras and Witt basis. In section three we will present a factoriza- 
tion for the operators {±idt — A) using an extension of the parabolic Dirac oper- 
ator introduced in [7] . For the particular case of the non-stationary Schrodinger 
operator we will present the corresponding Teodorescu and Cauchy-Bitsadze 
operators in analogy to |10| . Moreover, we will obtain some direct results about 
the decomposition of Lp-spaces and the resolution of the linear Schrodinger 
problem. 

In the last section we will present an algorithm to solve numerically the 
non-linear Schrodinger and we prove its convergence in L2-sense using Banach's 
fixed point theorem. 



2 Preliminaries 

We consider the m-dimensional vector space R"* endowed with an orthonormal 
basis {ei, ■ • ■ , e™}. 

We define the universal Clifford algebra C^o.m as the 2™-dimensional associa- 
tive algebra which preserves the multiplication rules eiCj + Cjei = —2Sij. A basis 
for C£o,m is given by eo = 1 and ba = ■ ■ ■ Ch,,, where A = {hi, . . . , hk} C 
M — {l,...,m}, for 1 < hi < ■ ■ ■ < hk < m. Each element x e C^o.m 
will be represented by a; = J^a^aga, xa € R, and each non-zero vector 
X = G a multiplicative inverse given by We denote 

by x'"^"''" the (Clifford) conjugate of the element x G C^o,mj where 

We introduce the complexified Clifford algebra C£m as the tensorial product 
C (g) Cio^m = ^w^Yl ^^^-4, e C, A c m| 

where the imaginary unit interacts with the basis elements via icj = eji,j = 
1, . . . , m. The conjugation in dm = C(X)C£o,m will be defined as w = lA^eA"^"'' 
Let us remark that for a, 6 G C£m we have \ab\ < 2'"|a||6|. 

We introduce the Dirac operator D = J^JLi ^j'^^i- It factorizes the m-di- 
mensional Laplacian, that is, — —A. A C4„-valued function defined on an 
open domain fi, u : f2 C M™ ^ C£m, is said to be left-monogenic if it satisfies 
Du = on n (resp. right-monogenic if it satisfies uD — on f2). 

A function m : f2 i— > C£m has a representation u — uasa with C-valued 
components ua- Properties such as continuity will be understood component- 
wisely. In the following we will use the short notation Lp(il), C'^(f2), etc.. 
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instead of Lp{^l,C£m), C''{ri,C£m). For more details on Clifford analysis, see 
0, [12], m or i. 

Taking into account [7| we will imbed R™ into R™+^. For that purpose we 
add two new basis elements f and satisfying 

f = f = 0, fft + ftf = 1, fe, + e,f = fe, + e,ft = 0, j = 1, ■ • ■ , m. 

This construction will allows us to use a suitable factorization of the time evo- 
lution operators where only partial derivatives are used. 



3 Factorization of time-evolution operators 

In this section we will study the forward/backward Schrodinger equations, 

{±idt ~ A)u{x,t) ^ 0, {x,t)€n, (1) 

where ft C M™ x IR+, m > 3, stands for an open domain in R™ x R+. We remark 
at this point that f2 is a time-variating domain and, therefore, not necessarily a 
cylindric domain. 

Taking account the ideas presented in ,1; and [7, we introduce the following 
definition 

Definition 3.1. For afunctionu G Wp{fl), 1 < p < +oo, we define the forward 
backward) parabolic Dirac operator 



D^,±,tu = {D + \dt±if)u, (2) 
where D stands for the (spatial) Dirac operator. 
It is obvious that Dx.±it ■ Wp{fl) Lp(ri). 

These operators factorize the correspondent time-evolution operator ([T]), 
that is 

iD,,±u)^u =. {±idt - A)u. (3) 
Moreover, we consider the generic Stokes' Theorem 
Theorem 3.2. For each u,v e Wp{fl), 1 < p < oo, it holds 



vdax.tu = / [{vDx-it)u + v{Dx,+itu)]dxdt 
n Jan 

where the surface element is d<Tx.t = [Dx + ^dt)\dxdt^ the contraction of the 
homogeneous operator associated to Dx^~it with the volume element. 

We now construct the fundamental solution for the time-evolution operator 
—A — idt- For that purpose, we consider the fundamental solution of the heat 
operator 

. ^ H{t) ( 
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where H{t) denotes the Heaviside-function. Let us remark that the previous 
fundamental solution verifies 

{~A + dt)e{x,t) = S{x)6{t). 

We apply to (|4]) the rotation t — > it. There we obtain 

(—A — idt)e{x, it) = —Ae{x, it) + dite{x, it) = 6{x)S{it) — -~i5{x)S{t), 

i.e., the fundamental solution for the Schrodinger operator —A — idt is 

e^(x,t) — ie(x,it) 

. H{t) f.\xn 

Then we have 

Definition 3.3. Given the fundamental solution e_ = e_(x,t) we have as 
fundamental solution i?_ = E^{x,t) for the parabolic Dirac operator D^^^it 
the function 

E^{x,t) = e^{x,t)D^^^it 



H{t) fi\x\^\ f ~X ef\W 



f +fM (6) 



(Anit)'^' '^^^ \ U J \2t ' ' 2t 

If we replace the function v by the fundamental solution in the generic 
Stoke's formula presented before, we have, for a function u G Wp{fl) and a point 
{xo, to) ^ 9r2, the Borel-Pompeiu formula, 

E_{x - XQ,t - to)da^^t t) 

on 



^ u{xo,to) + / E-{x ~ xo,t - to){Dx,+itu)dxdt. (7) 
Jn 

Moreover, if u e ker{Dx.+it) we obtain the Cauchy's integral formula 
E^{x - xa,t - ta)da.j,^t u{x, t) = u{xq, to). 



I an 

Based on expression ([7]) we define the Teodorescu and Cauchy-Bitsadze op- 



erators. 



Definition 3.4. For a function u 6 Lpi^l) we have 

(a) the Teodorescu operator 

T^u{xo,to) ~ / E^{x — xo,t — to)u{x,t)dxdt (8) 
Jn 

(b) the Cauchy-Bitsadze operator 

F_u{xo,to) = / E^{x - xo,t-to)dax.tu{x,t), (9) 
Jon 



for {xo,to) ^ on. 
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Using the previous operators, ([7]) can be rewritten as 

whenever v G Wp{Q), I < p < oo. 

Moreover, the Teodurescu operator is the right inverse of the parabohc Dirac 
operator D^^^n-^, that is, 

D^^^-uTu = / Dx-itE^{x ~ xo,t - to)u(x,t)dxdt 
Jn 



5{x — xo, t — tQ)u{x, t)dxdt 
u{xo,to), 



for aU {xq, to) e il- 
ia view of the previous definitions and relations, we obtain the foUowing 
results, in an analogous way as in 7j. 

Theorem 3.5. If v E Wp (dfl) then the trace of the operator F_ is 

ti{F.v) = Iv-^S^v, (10) 

where 

S-v{xo,to) = / E^{x - XQ,t - tQ)dax.tv(x,t) 
Jan 

is a generalization of the Hilbert transform. 

Also, the operator S- satisfies = I and, therefore, the operators 

2 2 2 2 

are projections into the Hardy spaces. 

Taking account the ideas presented in an immediate application is given 
by the decomposition of the Lp— space. 

Theorem 3.6. The space Lp{fl), for 1 < p < 2, allows the following decompo- 
sition 

Lp{n) = Lp{n) n ker{Dx,^it) © D^^^t (^Wp (ft) 
and we can define the following projectors 

P_ : Lp{n) ^ Lp{n)r\ker{D,^^,t) 

Q_ : Lp{n) [Wp {n) 
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Proof. Let us denote by (—A — idt)Q ^ the solution operator of the problem 

(—A — idt)u — f in n 

u — on dn 

As a first step we take a look at the intersection of the two subspaces 
DJ:,-^t (Wp i^)) and nker(D^,_,t). 



Consider u G Lp{Q) nkcr [D^.-it] HD^.-u yWp {^) j ■ It is immediate that 
D^^-itU = and also, because u G D^.-u I Wp {^) I, there exist a function 



o 1 

V ^Wp {^) with Dx^-itv = u and (—A — idt)v = 0. 

Since (—A — idt)Q^f is unique (see [TJ) we get v — and, consequently, 
u ~ 0, i. e., the intersection of this subspaces contains only the zero function. 
Therefore, our sum is a direct sum. 

Now let us u G Lp{n). Then we have 

-1 f °^ 

U2 = r'a;.-ji(-A - l9f)(, D^_,tU G ( Wp {^) 

Let us now apply D^^n to the function ui — u — U2- This results in 

= D^^-uu - D^^^uDx-iti-^ - idt)^'^Dx^_uu 
= D^^^itu - (-A - idt){-A - idt)Q^Da:-itU 
= Dx-itu - Dx,^ttu 
= 0, 

i.e., Dx^-itUi G Ver {Dx^~it)- Because u G Lp(fl) was arbitrary chosen our 
decomposition is a decomposition of the space Lp(fl). □ 

In a similar way we can obtain a decomposition of the Lp{Q) space in terms 
of the parabolic Dirac operator Dx^+u- Moreover, let us remark that the above 
decompositions are orthogonal in the case oi p = 2. 

Using the previous definitions we can also present an immediate application 
in the resolution of the linear Schrodinger problem with homogeneous boundary 
data. 

Theorem 3.7. Let f G Lp{n), 1 < p <2. The solution of the problem 

(—A — idt)u = f in VL 

u = on dfl 
is given by u = T-Q-T-f. 
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Proof. The proof of this theorem is based on the properties of the operator T_ 
and of the projector Because T_ is the right inverse of Dx-u, we get 

□ 



4 The Non-Linear Schrodinger Problem 

In this section we will construct an iterative method for the non-linear Schrodinger 
equation and we study is convergence. As usual, we consider the L2^norm 

I I/I I' = / [f7]odxdt, 

where [-Jo denotes the scalar part. 

Moreover, we also need the mixed Sobolev spaces W^'f'in). For this we 
introduce the convention 

^{x: {x,t) en} C M™ 
= {t: {x,t) en} C M+. 

Then, we say that 



u G w;^^'^{n) iff 



Under this conditions we will study the (generalized) non-linear Schrodinger 
problem: 

— AxU — idtu + \u\'^u — f in n (11) 
u — on dn, 

where = J^a l"-4p- We can rewrite (fTTj) as 

Dl_,,u + M{u) = 0, (12) 

where M{u) — \u\'^u — f. It is easy to see that 

u = -T^Q-T^{M{u)) (13) 

is a solution of ([T^ by means of direct application of to both sides of 

the equation. 

We remark that for u G W2''^{n)^ we get 

\\Dx,^.uu\\ = ||Q_T_Af(w)|| = ||r_Af(w)||. 
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We now prove that can be solved by the convergent iterative method 

Un = -T_Q_T_{M{u,,.i)). (14) 

For that purpose we need to estabhsh some norm estimations. InitiaUy, we 
have that 

||w„-u„_il| = ||T_g_T_[AfK_i)-AfK_2)]|| 

< Ci||Af(M„_i)-MK_2)||, (15) 

where Ci = ||r_Q_r_|| = ||r_||2. 

We now estimate the factor ||Af(u„_i) — M(w„_2)||- We get 

\\M{Un~l) - Af(u„_2)|| = |||w„-l|^li„_i - |u„-2|^Mri-2|| 
< |||m„_iP(u„_i - W„_2)|| + |||Wri-l - W„_2pW„_2|| 
< 2™+^||m„_i - W„-2|| (||w„-l|P + ||u„_2||||w«-l - Wn-2||) , 

We assume AC„ :— 2™+-^ (||w„_i|p + ||u,i-2||||w„-i — '«„_2||) so that 

\\Un - Un-l\\ < CllCn\\Un-l - Wn-2||- 

Moreover, we have additionally that 

||u„|| = ||T_Q_T_Af(u„_i)|| 

< 2™+iCi||7/„-i||3 + Ci||/|| (16) 

holds. 

In order to prove that indeed we have a contraction we need to study the 
auxiliary inequality 

2™+lCi||w„_l||3+Ci||/|| < \\Ur,-l\\, 

that is, 



+ ^ < 0- (17) 

The analysis of (fTT)) will be made considering two cases 
Case I: When ||w„_i|| > 1, we can establish the following inequality in 
relation to ^7}i 

||2 l^n-lll I ll/ll ^ n, ||3 , 11/11 



3 . 2™+i 2™+i ~ 2'"+-'-Ci 2™+-'-' 
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Then, from P7|) . we have 



\Un-l\ 



|Mn-l|| 



IK^ill , ll/ll 

1 



< 



6 • 2™+i 36 • 22™+2 2™+i 36 • 22™+2 

2 



< 



1 



\Un-l\ 



< 



1 



6 • 2™+i / ~ 36 • 22(™+i) 2™+i 



1 



1 



2™+i V 36 ■ 2™+i 



(18) 



If ll/ll < then 



1 



\U„-l\ 



6 ■ 2™+i 



where W 



ll/l 



In consequence, if 



W < ||m„_iII < 



1 



W 



6 . 2'"+! ■■ ^ ^ 6 • 2™+i 

then we have from p6p the desired inequality 

||Uri|| < ||w„_l||. 

Furthermore, we have now to study the remaining case. Assuming now that 
||m„-i|| < g.^i+i - W, we have 



|u„|| <2™+^C7i 



1 



6 • 2™+i 



< 



1 



6 • 2™+i 



W 



and ||m„_i|| < g^^i^i — W, ||u„_2|| < ^.^i+i — W so that it holds 



6-2'" + i 



\Un-l - Un-2\\ < 2 



1 



w 



Q . 2m+l 

With the previous relations we can estimate the value of JCn 

ICn = 2"+l(||li„_i||2 + ||7i„_2||||u„-l-U„-2||) 



< 2"+^ 



1 



< 3-2 



m+1 



6 • 2™+i 
1 

6 • 2™+i 



1 



6 • 2™+! 



- W 



W 



i_3.2™+iT4^< i, 
2 2 



(19) 
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which implies that 



\Un-2\\ < = 



1 



3 • 2™+i ■ 
Finally, we have that 

||u„-u„_i|| < /C„||w„_i - U„_2||, 

with fZn < \- 

Case II: When ||u„_i|| < 1, we can establish the following inequality 



\Un-l 



|4 ll'^n— 1| 



3 . 2^m+l 2"^+l 

Then, from (fT7|) . we have 



2"^+^C'i 2"'+-'- ' 



I ||4 l"r.-l||'' I ll/ll r, 

Pri-1 1 1 3.2m+l r 2^+1 ^ U 

^ (l|w„-l|P - 6.2"+0 - 3g . 22m+2 ^ 2"+l ' 

Again, if ||/|| < gg.^L+i then 



(20) 



1 



6 • 2™+! 



where W = 35:2^ " 5^ 
As a consequence, 



6-2 



6-2" + i 



<\\un-i\\ <J^:^ + W 



leads to ||m„|| < ||m„_i||. 

Again, considering now the case of ||u„_i|| < \J ^..^+1 — W , we obtain 

iKii <2™+ici(y=i;;^)'+ciii/ii <{^^, 

and 



W 



\Un-l - Un-2\ \ < ^J-^y^ - W 
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With the previous relations we can estimate the value of /C„ 
/C„ = 2"+i(||u„_i|P + ||m„_2||||u«-i-u«-2||) 



< 2"+^ 
= 3 ■ 2'"+! 



1 



6 • 2™+i 
1 

6 • 2™+i 



ly +2 



1 



6 • 2™+i 



ly 



i_3.2™+iw^<i 



which implies that 



|"n-2|| < R 



1 



3 • 2™+i 



Finally, we have that 



|u„-u„_i|| < /C„||w„_i - U„_2 



with /C„ < i. 



(21) 



The application of Banach's fixed point, to the previous conclusions, results 
in the following theorem 

Theorem 4.1. The problem hll]) has a unique solution u G ^^'^(rj) if f ^ 
L2{^) satisfies the condition 



< 



1 



36 • 2™+i ■ 



Moreover, our iteration method converges for each starting point 

o 1,1 

uq (^) such that 



koll < 



1 



6 • 2™+i 



with W 



1 NIL 

36-22(™+i) 2"»+i 
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